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Abstract. The solution of some equations involving functional derivatives is given as a series indexed by 
planar binary trees. The terms of the series are given by an explicit recursive formula. Some algebraic 
properties of these series are investigated. Several examples are treated in the case of quantum electrody- 
namics: the complete fermion and photon propagators, the two-body Green function, and the one-body 
Green function in the presence of an external source, the complete vacuum polarization, electron self-energy 
and irreducible vertex. 



PACS. 03.70. -f-k Theory of quantized fields ■ 
m Quantum electrodynamics 



11.15.Bt General properties of perturbation theory - 12.20.- 



1 Introduction 

Renornialization theory has been recently revitalized by 
the discovery a Hopf algebra that transforms the dreadful 
combinatorics of renornialization into a mechanical appli- 
cation of the Hopf algebra properties of rooted trees ||l|,||, 

^- n 

In the companion paper |^, Butcher's theory has been 
presented as an alternative way to describe the Hopf struc- 
ture of the algebra of renornialization. A particularly use- 
ful aspect of Butcher's approach is that solutions of non- 
linear differential equations can be written as a sum over 
rooted trees. 

In the present paper, Butcher's strategy is adapted to 
equations involving functional derivatives, that were first 
proposed by Schwinger ||] and will be called Schwinger 
equations in the rest of the paper. Schwinger equations 
are not commonly considered as a useful tool for compu- 
tation. The purpose of this article is to show that, by using 
series over planary binary trees, Schwinger equations can 
be turned into explicit calculation methods. 

The series we manipulate are indexed by planar bi- 
nary trees. So we first present some basic properties of 
planar binary trees. Then the solution of simple Schwinger 
equations are written as a sum over planar binary trees, 
with recursively defined coefficients. To make a compari- 
son with power series, the Schwinger equation would cor- 
respond to a differential equation for the sum of the series, 
whereas the formula we put forward corresponds to a re- 
cursive definition of the terms of the series: it does not 
contain so much information as the Schwinger equation, 
but it is more explicit to calculate the terms of the series. 

As an example, the full fermion and photon propaga- 
tors of quantum electrodynamics (QED) are written as a 



sum over planar binary trees. Renornialization is briefly 
described and applications are given for QED with an ex- 
ternal source, the vacuum polarization, the fermion self- 
energy and the irreducible vertex. 



2 Planar binary trees 

In contrast to Ref . H , we do not use rooted trees but pla- 
nar binary trees. Both can be drawn on a plane, but no 
permutation of vertices is allowed for planar trees. 



As an example. 



^and'\^ 



are two different planar 



trees, although they represent the same rooted tree. 

In planar trees, we distinguish two types of vertices: 
the leaves (which have no sons) and the remaining ver- 
tices (including the root), which we call internal vertices. 
In planar binary trees, internal vertices have exactly two 
sons. 

We now follow the notation of Loday and Ronco [g[. 
Planar binary trees have an odd number of vertices. We 
denote y„ the set of planar binary trees with 2n + 1 ver- 
tices. If t belongs to Ya, t has n -f- 1 leaves and n internal 



vertices. The number of elements of Yn is 



(2»)! 



the 



n\(n+iy.- 

Catalan numbers, which enter many combinatorial prob- 
lems and should probably be called Ming numbers ^ . 
If ti G Ym and ^2 ^Yn, the grafting of ii and ^2 is the tree 
t G Fm+n+i obtained by putting ii on the left of ^2 and by 
joining the roots of ii and ^2 to a new vertex that becomes 
the root of t. This operation is denoted by i = ti V ^2- For 



instance. 



grafting • and ^^ gives • V ^^ = \ X ■ 
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In the companion paper, rooted trees were graded by 
the number of their vertices. Here, planar binary trees 
have an odd number of vertices, and it is more natural 
to grade them differently: for each tree t, we define \t\ as 
the integer n such that t g y„. Thus, a tree t has 2\t\ + 1 
vertices. 

An essential property of planar binary trees is that 
each tree t different from • can be written in a unique 
way as ii V ^2, where ti and ^2 are called the branches 
of t. Moreover, grafting provides a recursive definition of 
planar binary trees |9|: 



r„= |Jnvr„_fc_i, Yo = {»}. 



(1) 



fe=0 



The notation Y^ V Yn-k-i means that all the trees of Y^ 
are grafted with all the trees of Yn-k-i- 

Planar binary trees have received much attention re- 
cently because of their relation to new algebraic structures 



3 Schwinger equations 

In this section, the solution of a linear Schwinger equation 
is given as a sum over planar binary trees. But we first 
introduce the concept of a functional derivative. 

A functional A{(l)) is defined loosely as a map sending 
a distribution to a complex number (see Ref.jll for de- 
tails). If i/) is a distribution, the functional derivative of 
A{cj)) in the direction t/j is defined as the limit for e ^ of 
{A{(t) + eij)) — A{(f)))/€. Finally, the functional derivative of 
A{4)) with respect to 4'{x), ^Jjll , is defined as the func- 
tional derivative of A[(j)) in the direction 5x, where 5x is 
the Dirac function 5x{y) — 5{y — x). 



3.1 Examples of functional derivatives 

A classical example is A[(f)) = j dxf{x)4>{x), giving easily 
j^rr^ — f{x). A further example, that will be useful in 
the sequel, is v4((/i) — J dxdyf{x,y)<j){x)<j){y). Then 

M(0) 



<5(/)(a;) 
S^A{(t>) 6 6A{(t>) 



5(j){x)5(j){y) 5(j){y) (50(x) 
More generally, if 



dyf{x,y)(j){y)+ j dyf{y,x)(j){y), 
f{x,y) + f{y,x). 



^(0) = j dxi... dxnfixi, . . . , Xn)(f>{xi) . . . 0(a;„), 



then 



(5"A(0) 



S<l){xi)-- ■54){x„) 



/ , f{Xa{l),.. .,Xo-(„)), 



where Sn is the set of permutations of n elements. 

In practice, A{(f) is often a Green function. Take the 
Green function defined by {Ax — 4>{x))A{(j)] x, y) — 6{x—y), 
which can be written A{(j)) = (Z\ — 0)~^. To calculate the 
functional derivative, we put A{(j) + eijj) = {A — (j) — eip)~^. 
The operator identity Y-^ = X-'^+X-\X-Y)Y-^ gives 
us A{(j) + eip) = A{(f)) + eA{(f))tljA{(j) + etp). Thus, taking 
the limit e — > 0, 



5A{(l>;x,y) 



dsA{(t); X, s)'ip{s)A{(f); s, y). 



If we choose now the distribution 'ijj{s) — S{s — z) we find 
6A{(l);x,y) 



6cb{z) 



A{(t);x,z)A{<j);z,y). 



(2) 



This identity will be used repeatedly in the sequel. 



3.2 A simple Schwinger equation 

As an introduction to the method of planar binary trees 
we consider the Schwinger equation 



X = A + F{X, 



sx 



dv{z) 



(3) 



where A is a functional of u, _F is linear in X and /^ ^ , 

' ov{z) ' 

and z is a variable over which F integrates. In an equation 
like (P), F is called the integral operator of the equation 
and A is called its initial data. 

A common example of such an equation is obtained 
when X = X{x,y), the initial data A{x,y) are the Green 
fun ction (A^ —v{x))A{x, y) = 5{x — y) discussed in section 
jj and the integral operator is 






dsdzX{x, s)f{s, z) 



SX{s,y) 
dv{z) 



for some function /(s, z). 

Such a Schwinger equation summarizes an infinity of 
equations that can be obtained by taking successive func- 
tional derivatives of Eq.(p|) with respect to v{z). 



SX , 

X = A + F{X,——) 



SX 



SA 



Sv{zi) Sv{zi) 

F{X, 



Sv{z) 
, 6X SX 



Sv{zi) Sv{z)' 
S^X 



Sv{zi)Sv{z) 



), 



aes^ 



When we take the n-th functional derivative of both sides 
of the (n-l)-th equation with respect to v{zn), the equa- 
tion gets an additional variable z„, and the chain rule is 
used to apply g_^^^ , to the right-hand side of the (ri-l)-th 
equation. 

If this is iterated to all values of n, we obtain an infinite 
system of nonlinear integral equations. This system seems 
difficult to solve because the n-th differential of X depends 
on the k-th differentials of AT for /c = to n + 1. 
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3.3 The series solution 

To write the solution of Eq. (|^) , we must introduce some 
notation. Sets of arguments will be often needed, so we 
write {z},j = z,, Zi+i, ■ ■ ■ ,Zj, {{z},j = \i j < i). Fur- 
thermore, if f{,{z}i,n) = f{zi, ■ ■ ■ ,Zn) is & function of n 
variables, then fi;{{z}i^n) is defined as the sum of n terms, 
where the first variable zi is shifted step by step from the 
first to the n-th position: 

fs{{z}l^n) = f{{z}l^n) + f{z2,Zl, {zjs^n) H h 

f{{z}2,n-l, Zl, Z„) + f{{z}2^n, Zl). 

For each planar binary tree t, we define an infinite- 
dimensional vector (j)(t), with components 0"(i), where n 
goes from to infinity. The n-th component is a function 
of n variables zi, . . . , z„. Now we define the initial data. 
For t — • , we take 0*^ ( • ) = A and 



Eq.(0), this defines (j>{t) recursively. The recursive defini- 
tion of (f>(t) is given explicitly by 

n 



fc=0 



•;2i 



5A 

Sv{zi)' 



(4) 



<Pl-^+\h; z,{z} k+i,n)), (8) 

for t ^ • and 0"( • ) is defined in Eq.(||). 

fn a quantum field interpretation, • represents the 
bare fields, V^ represents the interaction, and the sum 
over trees represents all the combinations of interaction 
that gives the full propagator. 

A proof of Eqs.(0) and (g) is given in the appendix. 



3.4 Enumeration 

If the initial data A is such that </<"(•) has only one 
term, as in Eq. dq) , the chain rule applied to the functional 
derivative gives a number of terms for 0"(i) that we de- 
note |0"(t)|. Equation (g) gives us the following recurrence 
relation for |(/>"(i)|: 



The most natural choice is to define 4>"( •) as l/n\ times 
the n-th functional derivative of A with respect to v(z), 
but this is not always the most economical choice in prac- 
tice. The only condition that we need for n > 1 is 



^•;Wi. 



Sv{Zn) 



0^(«;{4i 



(5) 



\r{t)\ = ^(n - k + f)|0'=(ii)||0"-^+i(t2)| 
Using the binomial identity 



When A(xj y) is a Green function of the kind discussed in 
section 3.f , it is not difficult to build such a 0( • ) from 
given initial data A. We can use Eq.ph to show that 

</)"(.) = A{:x,y), 
(/)^(«;zi) = A{x,zi)A{zi,y), 

</)"(. ;{z}i,„) = A{x,zi)A{zi,Z2)...A{zn,y), (6) 

satisfies Eq.(||) and condition (||). 

With this notation we can now write the solution of 
Eq.(pl) as 



fe=0 



a + k\ fb + n — k 



a+h+n+l 
a + b+1 



it can be shown that the solution of this equation is 



\4>-{t)\ = m 



2\t\ +n 
2\t\ 



X^Y.^'(t), 



(7) 



where ip{t) is an integer which depends only on the tree t 
(not on n) and is defined recursively by 

^(t) = (^(il)(2|i2| + f)<^(t2), 

ti and <2 arc the two branches of t. 



where t spans the set of planar binary trees. Moreover 
5X 



5v{zi) 



5''X 



5v{zi) ■ ■ ■5v{zn) 



Y,4>\t-zi), 

t 



o-e5„ t 



For each planar binary tree i, the vector 0(t) is calcu- 
lated as a function of the vectors (/)(ii) and 0(^2), where 
ii and t2 are the branches of t. Since </)(•) is defined from 



3.5 A compact notation 

To write a compact expression for the recursive definition 
of (/)(t) we define the deconcatenation of (zi, . . . , z„) by 



Z\(Z1,...,Z„) = ^(zi,...,Zj)(g) (Zj+i,...,Z„). 

ff z belongs to the vector space V , the map (f){t) acts on the 

00 
tensor module (Fock space) T{V) — ^V®^. Sometimes, 

n=0 
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as in QED, (f) is defined on T{V) x M, wfiere M is a 
fixed vector space, for instance M = V^ for the photon 
propagator. 

We define the operator d{z) by 

d{z)<p"{t;zi,.. .,z„) = (/)J.+^(i;z,zi,.. .,z„). 

The recursive definition of (j) becomes 

0(i) =Fo{ld(g> d{z)) o {(j){ti) ® (j){t2)) o A, 

where F{a ® &) = F{a, b). 

It would be interesting to find a family of equations 
such that, for any (j) from T{V) to C, there is a member 
of the family of which </> is a solution. This would gener- 
alize Butcher's density theorem[0], and would provide a 
general class of equations that would be satisfied by the 
renormalized Green functions. 



3.6 Algebra structure 

In the case of rooted trees. Butcher [|4[ has defined a 
group structure of Runge-Kutta methods that Hairer and 
Wanner [^ interpreted as a composition of Butcher se- 
ries. A similar approach can be used for planar binary 
trees. A powerful aspect of Butcher's approach is that al- 
gebraic operations are defined on two spaces at the same 
time: the space of Runge-Kutta methods, and the space 
of maps over trees. The same strategy will be used here, 
and the operations will be defined on the space of integral 
operators and on the space of maps over planar binary 
trees. 

We start with the addition. If we have two Schwinger 
equations X = A + F{X, ^) andY = B + G{Y, ^), the 
addition of the integral operators is H = F + G and the 
addition of the maps 4> (corresponding to the first equa- 
tion) and tp (for the second equation) is defined recursively 
by x"( •) = </)"(•) -h -0"( • ), for the initial data A + B 
and 

X"(i;{z}i,„) = 

n 

J2 F{x\h- {z}i,fc), xr'+'fe; z, {z},+i^„)) + 

n 

J2 G{xHti; {zKk),xT'+\h- z, {z}fe+i,„)). 

fe=0 

This addition defines clearly a commutative group struc- 
ture for the integral operators. It defines also a commuta- 
tive group structure for the space of maps, where the unit 
element is (/)(t) = for all t, and the opposite of (j){t) is 

^(i) = -(-l)l*l</'(t). 

Multiplication by a scalar A is defined similarly. An 
integral operator F becomes XF, and the corresponding 
map (j){t) becomes Al*l0(i). Notice that maps are not equiv- 
alent to integral operators since they contain also the ini- 
tial data. The present definition of the multiplication by 
a scalar corresponds to the case where the initial data are 



not changed. If the initial data are also multiplied by A, 
then (j){t) becomes AI*I+V(^)- 

This addition is useful when a Schwinger equation is 
the sum of various terms. Addition and multiplication by 
a scalar is probably all we need to renormalize quantum 
electrodynamics. However, we can proceed with Butcher's 
strategy and define another operation coming from a com- 
position of solutions. If we start from two Schwinger equa- 
tions X = A + F{X, ^) and Y = B + G{Y, ^), the com- 
position of the solutions is defined as the Y obtained with 
the initial data B = X. 

It is shown in the appendix that, if x is the map corre- 
sponding to Y (i.e. F = Et x(0), then x(i) = </>(*) +V'(i), 
where 0(i) is the map associated to the equation for X 
and V'" (t) is given by ■0" ( • ) =0 and 

n 

r{t; {z}i,n) = Y. G(0^ii; {2}i,fc) + i>Hh; {z}i,k), 

k=0 

rr-'^Hh; z; {z}k+i,n) + ^r'+'(*2; z; {z}fe+i,„)). (9) 

This defines a product of integral operators and of maps. 
In the first proof of Eq.(P) given in the appendix, the 
integral operator corresponding to this product is con- 
structed. Notice that this operator acts on vectors (x+y) ■ 
This product has a unit element (given by G = 0). 

In the next section, the present method will be applied 
to the example of quantum electrodynamics (QED). 



4 The case of QED 

We work in the flat Minkowski space with a diagonal met- 
ric g (the diagonal is (1, —1, —1, —1)). The electron charge 
is e = — |e|. Repeated indices are summed over. 

In 1951, Schwinger B devised coupled equations in- 
volving functional derivatives of S{x, y; J), the full fermion 
propagator of QED in the presence of an external electro- 
magnetic source J^{x): 

[ag^, - (1 - ^)d^d,]A''ix; J) = -J^{x) - 
ieiY['y^^S{x,x;J)], 



^I'^d^ 



e-f'^A^{x;J) -f ie7^ 



5J^{x) 



S{x,y;J) = 
Six~y).{lO) 



Building on a work by Polivanov |l6| , Bogoliubov and 
Shirkov JTtI transformed this equation into a Schwinger 
equation coupling the full fermion propagator S{x, y) with 
the full photon propagator D^^{x^y): 

[Dg^- - (1 - i)d^d'']D,p{x, y) = g>;5{x - y) - 

'-^Ar^-^^^]D.p{z.y;Al (11) 



(12) 



5A,{z) 

[il^d^ - TO - e7^v4^(2;)] S{x, y; A) = S{x - y) + 

SS{x,y;A) 



ie / d'^zj''D^p{x,z;A)- 



SApiz) 
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where A{x) is now an external electromagnetic field. As branch is on the left and the black (fermion) branch is on 
explained in Ref. |0[, Eqs.(ll) and ( |l^ ) are not completely the right. 



equivalent to Eq.(M), they are valid in the limit A — For example, the fermion trees are 

(J = 0), which is the standard case of QED. 

Multiplying E g ([ill ) by the bare photon propagator, o . V V V Y 

D^,,^{x,y) and Eq.(|12D by the bare fermion propagator in « \j \j \/ 

the presence of A, S°{x, y; A) = [i"f^d^ — m — e7^A^]~^, 

we obtain our starting Schwinger equations: and the photon trees are 

D,,{x,y;A)^D%{x,y)~iejd^zd^z'Dl^{x,z) o Y \f Y/ XX' '^ 

1 ~T7~"T3\~J ^^'v{z , y; A), Frabetti has calculated the number of trees with 2p + 

„ 2q+ 1 vertices, p+ 1 white leaves and q + I black leaves 

S{x, y; A) = S^{x, y: A) + ie / d^zd'^z'S°{x, z: A) as 



An . _, 6S(z,y;A) _ {p + q)l (p + g+1)! 

In principle, these equations fully determine S{x, y; A) and The following recursive definition gives (/>(i) in terms 

D (x y A) of '/'(^i) ^i^d (pit'), where t^ and t* are the branches of t: 



4.1 The tree solution 



The method of the previous section is now used to write 



n n — k „ 

0"(i';x,y;{A,z}i,„)=i^^ / dSdV 

J n J. / n "^ 



in—k—k I 



" {•■,x,z;{X,z}i^n-k-k') 



the solution of Eq. (|T^) . Since S'"(x,z;A) depends on the J^fpxx' {ti; z, z' ; {X, z}n-k-k'+i.n-k') 

external potential A(a;), the small extension presented in ,fe'_i_i/ ,_ .\' ' r\ i ^ n 

the appendix is required. All quantities will be taken at 9s \ 2^ ' U^ ^ ^i ^ /n-fc'+i,nj, v 

A(a;) = 0, so the external potential will not be mentioned ,r, / o V^ /" j4 .4 / n / 

(e.g. Six, y) means Six, y; A) at A^ 0). "^/^-(^ ' ^' V'^ i^' ^^i.") = -^2^J d zd z D^^ix, z) 



fc=0' 



The notation {A, z}ij = Xi, Zi, Xi+i, Zi+i, . . . ,Xj, Zj en 
ables us to write the solution as "^Avl^i! ^ i J/! l-^i ^}i,fc) 



tr [7 Vr '+' (i* ; ^, ^; A'^', {A, z}fc+i,„)] . (17) 

This recursive definition is completed by giving the 
5Six,y) _ Y^ 2|i|+ii,i/,.. .\ N components of 0"(») and (/)"(o): 



t 

A4 ^ ^ - Ve2|*l+V'(i-;x,2/;AiZi), 
oAx.izi) ^ 



£'M''(a;,y) =5I^^'*''?^"''(*°'^'y)' (15) 0;^^(o;x,y) = for n>l, 

a"S'('t w^ ^^ ^^ ,, 'P^i»;x,y) = S^i^,v), 

" '^K^^y) ^ V^ \^p2|t|+n n A n 

6Ax,izi)---SAx^iz„) ^^ (t> i»;x,y;XiZi) = S ix,zi)-f 'S izi,y), 

"^%r*f--.„-n .; ^ ri';x,y;{X,z}i^n) = S"ix,zi)-f^^S°izi,Z2h^---- 

(•ji T) c \ T *^ \^Tii y)' 

Mai (2:1) ■ • • SAx„ (z„) J^ ^ In practice, the double sum of Eq.(|l|) is replaced by: 

<Pl^it°;x, y; {A, z}<.(i),^(„)). 0"(i'; a;, y; {A, z}i,„) = 

Another change in the notation is that trees have now two '-' i^,z\)^ ^cp (i ; zi, y; {A, z|2,n) + 

colors. The map (/)(i) has a supersymmetric flavor because "^ /■ 

it has a fermion component for the full fermion propagator *Z^ d zd z S (x, 2)7 (/)_;^;i,, (i^; 2^, 2: ; {A, zjifc) 

and a photon component for the full photon propagator. It k=o 

is convenient to transfer the fermion/photon index from ip d/ir'^^^itl- z y; X'z' \X z}k+i „)■ (18) 

to the tree t. Thus, 0(t°) is the photon component of (pit) 

and (/)(t*) is its fermion component. Futhermore, 0(t°) will This is still a recursive definition, but it uses now the 

be drawn with a white root and (pit*) with a black root, smaller components of cp for the same tree t*. Notice that 



All trees are built as t^ V t*, where the white (photon) the first term of Eq.(18) is absent if ri = 
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Three remarks can be useful at this point. Firstly, con- 
sidering all quantities at A = 0, we find as in section \i.l\ 
that SS"{x,y)/SAx{z) = eS°{x,z)-f^S°{z,y). Therefore, 
in the definition of ^" ( • ) , a factor e" was suppressed and 
transferred to the solution Eqs.(O) and (p^). It must be 
checked that this is compatible wit h re normalization. Sec- 
ondly, it will be shown in section 6.1 that, from Furry's 
theorem, the components (jf^(t°) are zero when n is odd. 
This reduces the sums in Eqs.(16), ( [l7| ) and ( |l8|) to the 
even co mp onents of 4>{t1). Finally, the vector space V of 
section 3.5 has now become the space {0, 1, 2, 3} x R"'. 



The next term is 



(Y) 






4.2 Diagrammatic interpretation 

The recursive solution of the previous sections can be il- 
lustrated in the usual diagrammatic language: 



/ 



< o o o 


V ; ; 



JJJ^ 



\ 



\ 



I 



In other words, 0"( • ) is a propagator with n dangling 
photon lines {n starts at 0). More generally, all the com- 
ponents (i>^{t°) or 4)^{t') has n dangling photon lines. To 
see the action of the recursive equation for the photon, 
observe that, in Eq.([l7|) the fermion extremities of t\ are 
closed on an additional photon line by a bare vertex (on 
the left), and each dangling photon line of t* is linked in 
turn to the photon extremities of f^ (on the right). Dia- 



^rammatically: -^^ ^^~'^-^~" 



The next term is 



(Y) 






v 



/ 



For the fermion propagator, in Eqs.([l6D and (|18|), each 
dangling photon line of t* is linked in turn to the right 
photon extremity of f^ ■ The left photon extremity of i2 is 
linked to the left extremity of the electron propagator t\ 
by an additional bare vertex. 

If we neglect the first term in Eq.(p^ we can write 



diagrammatically: 




4.3 Diagram enumeration 

As a warm-up exercice, we can calculate the number of 
diagrams in the component (j)'^{t), that we denote |0"(i)|. 
From Eqs.(|l6|) and ( p7[ ) we find the equations for |0"(i°)l 

and \(f>'"{t')\ 

n n — k 



fc=Ofc'=0 



w'^{n\^Y.^n-k + i)\4>\ti)u 



ri-fc+l 



(^2)1, 



fc=0 



With0"(.) = l, 0"(o)=J„,O; 

The solution of this recursive equation is 



|0"(i-)|=^(t) 



l0"(i°)l=^W 



2|t| 



2\t\ +n- 1 



where ^p{t) does not depend on the colour of t and was 
defined in section 3.4. 



This can be used to calculate the total number of dia- 
grams for the electron or photon propagators (the number 
is the same) contributing at e^". Let us define 



teY^ 



teY„ 



Using Eq.(P for n > we find 



n-l 

Sn = EE'^(^l) E (2|i2|+l)^(i2) 

fc=0|ti|=fc \t2\=n-k-l 

n-1 



^(2fc + l)sfcs„_fc_i. 



fc=0 

The starting value is sq — 1. For n — 0,1,2,3,4,5 this 
gives us 1,1,4,27,248,2830, in agreement with Refs.(l|, 
n3. The generating function y{x) for the sequence s„ sat- 
isfies the differential equation 2x^yy' + xy^ — y + \ — Q 
with y(0) = 1. 

This enumation takes into account neither symmetry 
nor Furry's theorem, which says that |(?!)"(i°)| = if n is 
odd. The main point of this enumeration is to show that 
each tree represents the sum of a large number of diagrams 
when |i| is large. This may prove useful for practical cal- 
culations. 
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4.4 Fourier transform 

In applications, it is often convenient to work in the k- 
space. The Fourier transform of (/i(i; x,?/; {A, z}i^„) is de- 
fined by 

V'(i;g,<7';{A,p}i,„) = / d'^xd'^yd'^zi . . .d'^Zn 

This corresponds to outgoing momenta pi along the dan- 
gling photon lines. If this is introduced into Eqs.(|lq) and 
(y), we find 

4>{t;q;{\p}i,n)- 

The full fermion and photon propagators in Fourier 
space are 






Here (j){t) satisfies the recursive relation 

^"(i-;q;{A,p}i,„) = 5°(g)7^^ 

r-\t';q+p,;{X,p}2^^) + 



^E 



^ P cO/„N. A 7 



S^\qh'^tx'{tl;p;{X,p}i,k) 



(27r)4 
^;:,(r;q;{A,p}i,„) = -* E / TCT^maC?) 



fe=0 






fc=0 



{2^Y 



tr[7^0^-^-+i(i';p; A', -q - P,, {X,p}k+un)] , (20) 

where we have noted Pk = pi + - ■ ■ +Pk and with the initial 
data 



4>U'''q) = DU'i) 



gtiy - (1 - l/0'?^9i'/9^ 



q'^ + le 

0M-(°;9)=O for n>l, 
~^\.-q)^S''{q)^{rq^.-m + ier\ 
~4>\.-q-Xi,p^) ^ S\q)^^^S\q+p,), 
I'-i.-q- {A,p}i,„) = 5"(q)7^i^0(g -I- p^)-1^- • • • 7^" 

S°{q + pi + ---+Pn). 



As for the real space case, the first term of Eq.(|l9|) is 
absent for n = and 

0'i;(t;g;{A,p}i,„) = <^"(i;g;{A,p}i,„)-f 

<^"(t; g; X2,P2, Ai,pi, {X,p}z,n) H h 

'^"(^;9;{A,p}2.n,Ai,pi). 

Again, Furry's theorem enables us to restrict the sum to 
the even components of 0p^(t°). 



5 Renormalization 

In Ref. Ql, it has been shown that the renormalization 
through Hopf algebra created by Kreimer could be done 
also within the Butcher approach using a difference of 
Runge-Kutta methods (bare quantities minus countert- 
erms). The same method can be tried here. 

li X = A + F{X, gj(^\ ) is the equation to renormalize, 
and if R represents a renormalization prescription (see 
Rcf . H ) , then the renormalized field Xr is given by 



^ -^M^J, 



t 



where (t)R{t) is defined recursively from the integral oper- 
ator i^ - i? o F by 0^( • ) = 0"( • ) and 

n 

rR{t;{zh,„) = Y^ 

fe=0 

i?[F(4(ti;{z}i,fc),(/)S^'^-+i(i2;z,{z}fc+i,„)) . 

Within this approach, renormalization looks rather sim- 
ple. At each stage of a calculation, use the renormalized (j) 
instead of the bare ones, and subtract the singularity. This 
rule seems promissing for the case of QED, where </)"(•) is 
regular and needs not be renormalized, and where the re- 
cursive operations described by Eqs.(|l9|) and ( po| ) involve 
only a single loop (a fermion loop for the photon propaga- 
tor and a photon loop for the fermion propagator). Since 
0|.(ii) and (/)^~'^^^(i2) are regular, there is no subdiver- 
gence, and thus no overlapping divergence problem. 

Following Kreimer , the renormalization group equa- 
tions are a consequence of the fact that F — R{F) — 
R'{F)-R{F)+F-R'{F) is independent of the parameters 
oiR'. 

Clearly, this presentation of renormalization is far too 
sketchy and further investigation is required to include 
renormalization prescriptions and calculation of countert- 
erms, and to check that the results are the same as in the 
standard theory. 



6 Applications 

In this section, the previous results are applied to Furry's 
theorem and the Ward-Takahashi identities. Finally, a sum 
over trees is defined for the two-particle Green function. 



6.1 Furry's theorem 

Within our approach, Furry's theorem implies </>" {t° ) ~ 
for odd n. To show this, remark that, in Eq.(EOh, the 
integral 



d^p 

(2^ 



tr[7^0J-'=+i(t-;p; A', -q - Pk, {A, p} ,.+!,„] 
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is a fermion loop with n — k -\- 2 external photon lines. 
According to Furry's theorem, this loop is zero when n — k 
is odd. From its explicit definition, </)"( o ) = is zero if n 
is odd (in fact, it is zero if n > 1) . We reason recursively on 
the number of vertices of i° . If (^" (i° ) = for odd n and t° 
with up to 2 A''— 1 vertices, take t° with 2N +1 vertices. In 
Eq.(EO), the integral over p is zero if n — fc is odd (Furry's 
theorem) and (t>*'{ti) is zero if k is odd (because t^ has less 
vertices than t). Thus, for the left-hand side of Eg. (po[) to 
be eventually non-zero, k and n — k must be even, so n 
must be even. Therefore, (t>^{t°) = for odd n. 



We define di by 

9i^5;„ (t; q; {A, z}i,,_i, Ai, k, {A, z}i+i^„) = 

3 

X! ^M^'l"' (*; 9; {^^ z}i,,^i,^i, k, {A, z}i+i^„) == 

^^7.-1 (i; g; {A, ^li.i-i, {A, zj^+i^n) - 

^'^7-1 (<; 9 + ^; {A, z}i,j-i, {A, z}i+i,„), 

where the last line is the Ward-Takahashi identity. If we 
write the usual expression for the boundary operator 



6.2 Ward-Takahashi identities 

Within the present approach, the Ward-Takahashi identi- 
ties take the form 

Y,k^4>l+\t';q;fi,k,{X,ph,,) ^ rit';qA\p}i,n) - 

rit';q + k;{X,p},,„). (21) 

Again, these identities will be proved recursively over 
the number of vertices of t*. An elementary calculation 
shows that Eq.(|l|) is true for f = • (Ref.|§, p.238). 
Using Eg. (p9|) , we can write Eg.(pi|) as 

E <50"(i*;g;{A,p}i,„) ~ 
^M XA n„\ — — = <?^ (* ;'7;{A,p}i,„) 



6'i;„ (t; q; {A, z}i,„) = ^(~l)^9,^'^„(t; 9; {A, z}i,„), 



1=1 



we find by explicit calculation that d^ = 0. This was 
proved for the fermion component, but d^ = for the 
photon component too, because of Furry's theorem. No- 
tice that, for n = we have put d(ffi{t) = 0. 



6.4 Two-particle Green function 

In Ref.pj Schwinger proceeds by giving an equation for 
the full two-particle fermion Green function: 



SA^ik) 



r{t';q + k;{X,p},,n). (22) 



The last eguation can be proved by applying the differen- 
tial operator ^ ^ij- sa (k) ^^ both sides of Eg.dl9). Con- 
sider first the second term of Eg. dl3) . The differential op- 
erator acting on 0^^, (tj) adds a photon fine to (j)\y{t°), 
which gives zero by Furry's theorem. The differential op- 
erator acting on (f/^ '^^{t') will transform according to 
Eq.(p2|), which is valid because t* has less vertices than 
t* . Carrying out the sum over k, we see that Eq.( P2[) i s 
true for the second term of the right-hand side of Eg.(|19D. 
Consider now the first term of Eq.(|l9|). For n = it 
does not exist, so Eq.(|2|) is proved for n = because it 
is true for the second term. The recursion runs now over 
n, the functional derivative operates either on S^{q) or on 
4)'^-^{f]q + pi;{X,p}2,n)- Since Eq.(||) is true for both, 
and for the second term, it is true for (/'"(t*; q; {A,p}i_„). 



6.3 Homology and Ward-Takahashi identities 

The previous Ward-Takashi indentities induce a boundary 
operator on 4>(t). In that case, we need to start from the 
fully symmetrized form of 0: 

creSn 



L 1 SJfj^{xi). 

S{xi,x2;yi,y2;J) = S{xi -yi)S{x2,V2\J) - 
5{xi -y2)S{x2,yi;J). 

From the identity [|l^ 

we go from the J to the A variable (in the limit A=0) 

il^d^i-^ - m - e-i^" A^,{xi)]S{xi,X2;yi,y2; A) = 

5(xi - yi)S{x2,y2; A) - S{xi - y2)S{x2,yi;A) + 

■ f ,4 uSS{xi,x2;yi,y2;A) 

le j d Z7^ JaTz) G\t^{z,xi). 

It remains to multiply by S^ and to take A — Q. This 
yields 

S{xi,X2;yi,y2) = S°{xi,yi)S{x2,y2) ~ S°{xi,y2)S{x2,yi) 



+ie / d''zd^zS"{xi,z)-/'' 



6A^{z') 



-Gx^{z',z). (23) 



In Eq.(p3|), S{x, y) plays the role of initial data, whereas 
it is the solution of Eq.(n3). Therefore, we are making a 
composition of solutions, which was considered in section 
^.6| . The situation is not exactly the same as that of sec- 
tion 3^, but the proof is similar (only more cumbersome) 
and we obtain an expression for the two-particule Green 
function as a sum over planar binary trees 

S{xi,X2]yi,y2) = ^e^l*lx°(^;a;i,a;2;2/i,2/2)- 



Christian Brouder: On the trees of quantum fields 



According to the rule of composition, x is the sum of two 
terms: 

S°{xi,yi)(j)°{t';x2,y2; {A, z}i,n) - 
S°{xi,y2)(l)°{t';x2,yi; {A, z}i,„) + 
il;°{t;xi,X2;yi,y2;{>^,z}i,n), 

and ip itself is given by 

■0"(i; a;i,X2; 2/1,2/2; {A, z}i,„) = S°{xi,zi)j^^ 
x"~^it;zi,x2;yi,y2;{\z}2,n) + 



subtree of t that has Si as a root, and Ui the subtree of t 
that has Si as a leaf. For example: 



^m-vv. 



This operation will help calculating the vacuum polar- 
ization and the self-energy. If we define Y by 



n ^ 

J2 / d^zd^z'S'>{xi,z)j^cf>'^^{tl;z',z;{X,z} 

1 n '^ 



l,k) 



fc=0 
n-k+l 



Xs ~^ {h; z,x2;yi,y2; >>, z\{X, z}k+i^n)- 



Equation ( |24D is not very simple, but it provides a 
way to calculate recursively all orders of the perturbation 
expansion for the two-particle Green function of QED. In 
that sense, it is not so complicated. 



7 Vacuum polarization, self-energy 

To calculate vacuum polarization and self-energy, we in- 
troduce a further operation on planar binary trees. Let 



<t) 



P{t) =^U,®Vi 



(24) 



where n{t) is an integer recursively defined by 

n(.)=0, 
n{t) =0 if t == ti V • , 
n{t) = l + n{t2) if t^ti\Jt2, t2^», 

and the planar binary trees Ui and Vi are determined by 

P(.) = 0, 
P(t) =0 if t = ti V . , 

n{t2) 
P{t) ^{hV •)®t2+^{tl\/ Ui) (g) Vi 

4 = 1 

if t = hVt2, t2^ •■ (25) 



The trees Ui and Vi in Eq.(p5[) are generated by Eq.(p4[) 
for t = t2- 

As a more graphical definition, for a tree t, we consider 
the path starting from the root and climbing up the tree 
by taking, at each vertex, the right branch. This path 
terminates at the extreme right leaf of the tree and goes 
through n{t)-\-2 vertices (including the root and the leaf). 
For each vertex Si along that path, excluding the root and 
the leaf, we cut t into two trees Ui and Vi, where Vi is the 



it is proved in the appendix that Y — ^^ 4''^{t), where 

V'''(.) = o, 

1 1 "'*^ 

V'°(t) = ^ij(7y(/W^ij(7y -E^°(-0^°(«.))(26) 

In this equation, m^, Vi and n{t) are determined from t by 
Eq.(pl). 



7.1 Vacuum polarization 

From the general formula, we deduce that the vacuum 
polarization is given, in terms of the map for the full 
photon propagator, by 

t 

with ■(/;0( o ; g) = and 
^l^t°\ct) = (7"9a - m)0O^(r ; g)(7'3g^ - m) - 

n{t) 

(7-gc. - m) ^ 0tK; '?)ff''''>°>«; «)• 
j=i 

Incidentally, it can be noted that the full vacuum po- 
larization n^^ (q) is related to the fermion loop by (Ref. |gl[ 
p.477): 

n'^iq) ^-^Y: e^l'l^^ / ^M^'4>\t';p; ^^, -q)] . 



7.2 Self-energy 

Similarly, the self-energy is given, in terms of the map 
for the full fermion propagator, by 

r(g)=^e2|*IV'°(i*;(?), 
t 

with ?/'0( • ; g) = and 

^°it';q) = il'^qa - m)<j,°{t';q){j^q0 - m) - 

nit) 

{^^q^-m)Y,<lP{u:;q)^\v:-q). 
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From the formula for the self-energy, we can deduce 8 Interaction with an external field 

the complete one-particle irreducible three-point function 

r^{p,p') [pi[, p. 335. We reintroduce the external potential In this section we come back to the original Schwinger 



to write 

S{x, y; A) = ^(i7^9,a - m - e-f'^A^^ix)) 
t 

(t>{t*;x, y- A){-i^^'^y^. - m - eYA^^iy)) - 

^(z7^5^A -m-e-i^Ax{x)) 
t 

J2 J dz'^{ul-x, z'; A)V'«; z', 2/; A), 

i—l 

where d yi^ acts on the left. In the real space, r'^{x, y; z) 
is given 

rvi \ SS{x,y;A) . . „ 

1 (x,y;z) = — :r-. — r^ — for A = 0. 



SA^iz) 



Therefore, 



t 
with 

■il;^{t';x, y; v, z) = -^^{z - x) 

<iP{t'-x,y){~i^^'"dy.-m)- 

{i^^d.,.-m)4P{t';x,y)Y5{z-y) + 

(17 d^\ — ■m)(j)^{t'; X, y; v, z)(— 17*^ d yi^ — m) + 

"(t) „ 
j-'Siz -x)J2 dz'cfPiu^x, z')4,\v:-z\ y) - 

i=i •' 
"(*) „ 
{ij^d^x -m)^ dz'(t)°{u*;x,z')ip^{v*;z\y;i^,z) - 

ih^d^x -m)^ / dz'(l)^{u']x,z';iy,z)il;°{v';z',y). 

i—l 

In Fourier space, this gives us 

t 
with 

ij\t';q;u,p) = -7>°(i*; «)(*7^(9p +P^)-m) - 
{l-f^qx~m)4P{t'-q)Y + 
{il^qx - m)(j?-{t';q]v,p){i-i^'{q^, + p^,) - m) + 

n(t) 

j''J2^%u:;q+p)^p'{v:;q+p)- 

i=l 

n{t) 
{il^qx -m)'^(j)°{u';q)ip^{v*;q;v,p) - 

i=l 
"(t) 

(«7^<7A -m)^ </>^«; <7; v,p)^J°{v*\q + p)- 



equation, because the presence of an external source is 
a convenient way to represent the nuclei in the QED of 
matter. 

Starting from Eq.(|lO|), we multiply by the correspond- 
ing bare Green functions and we introduce A'^{x; J) into 
the second equation to obtain 



S{x, y- J) = 5"(x, y)-ej d^zd^z'S\x, z)-i^'D%{z, z') 

r{z')S{z, y; J) - le^ f d'^ zd* z' S° (x , 2)7'' 

D%{z,z')tr[j''S{z',z'-J)]S{z,y;J)- 
,SS{z,y;J) 



le f d'^zS%x,z)J''- 



SJ^'{z) 



(27) 



Here, the Schwinger equation is a sum of three terms. 
The first term is simply the classical interaction with the 
external source J'^{z2), it can be solved by defining a bare 
propagator in the presence of this source 

S°ix,y-J)-^ = i7^5,, - m + 67^ f d'^zDl^{x,z)r{z). 

Equation ( pTf ) becomes now 

S{x, y; J) = S"ix, y; J) - ie^ [ d^zd^z'S"{x, z; J)7^ 



Dl,{z,z')tT[rSiz',z';J)]Siz,y;J)- 
le d zS {x,z;J)-f^'—^——. (28) 

This equation is solved by the usual methods, and the 
recursive definition of (/>(i) is 



r( • ;x,y; {A, z}i,„) = (-!)" y d*si . . .^^^^^"(x, si; J) 
7^"i?°„,Js„,z„)5°(s„,y;J) 

n n—k „ 

0"(i; X, y- {A, 2}i,„) = -^Y.Y. '^'^^'^' 

fc=Ofc'=0'^ 

</."-'=-'='(. ; X, z- {A, z}i,„_fe_fe07^^".(^, A 
tr[7'^0''(ii; z', z'; {A, z}„_fc_fc/+i,„_fc')] 
0'' (^2;2:,?/;{A, z}„_fc/+i,„) - 



;^ / (i''z(?!)''(«;x, z;{A, z}i,fc) 

fc=0 



iA'e ''^H^2;2^,2;;M,^,{A,z}fc+i^„), 

where the last term is non zero only if t has the special 
shape i = • V ^2 . 

According to Schwinger y| , the full photon Green func- 
tion in the presence of an external current J is given by 
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the functional derivative of A{x; J) with respect to J{y). 
Therefore 



Dx^,{x,y;J) 



&Ax{x-J) 
SJt^iy) 
DUx,y) + 



ie / dzZ?^^(a;, z)tr [7 



^SS{z,z;J)-. 



SJt'iy) J 
ie^^ / dzDl^{x, z)tr[j''(l)'^{t; z, z; n,i 



It is also possible to start directly from Eq. ( p7| ) and to 
write a tree solution of this equation using bare fermion 
Green functions. Here, the strong field case was treated 
because it is probably more interesting for applications to 
solid-state physics. 



9 Conclusion 

A method was presented to write the solution of some 
Schwinger equations as a series over planar binary trees. 
In quantum field theory, it is common to expand over the 
number of loops or to use integral equations relating, for 
instance, the full propagator to the full vertex. The first 
method gives explicit results but becomes very complex, 
and hundreds of diagrams must be built and calculated af- 
ter the first few terms of the perturbation expansion. The 
second method is formally powerful but not very explicit 
because an n-body Green function is expressed in terms 
of an unknown (n-l-l)-body Green function. The present 
approach is a way to mix these two methods to obtain 
an explicit recursive formula for the propagators and its 
functional derivatives. 

The main point of the method is that explicit recur- 
sive expressions can be given for the solution of Schwinger 
equations. Because of the recursive structure, the results 
obtained at each step can be reused for the next step. 

The present paper is only a first exploration of the 
method of series indexed by planar binary trees, and much 
work remains to be done to explore its algebraic properties 
and its applications. 

Two kinds of applications were presented here. On the 
one hand, a series indexed by planar binary trees was given 
for various physical quantities (full fermion and photon 
propagators, full two-body Green function), and a method 
was given to deduce from this a series for vacuum polariza- 
tion, fermion self-energy and irreducible vertex function. 
Although the formulas may be a bit cumbersome, they 
are derived and proved easily. On the other hand, the re- 
cursive nature of the terms of the series is well suited to 
prove properties to all orders of perturbation theory. This 
was done in the case of the Ward-Takahashi identities. 

The present work can be expanded in various direc- 
tions. Other field theories can be investigated, as well as 
many-particle Green functions. The present treatment was 
restricted to classical electromagnetic sources, it could be 



worthwile to study the case of anticommuting fermion 
sources. Furthermore, planar binary trees could be used 
to solve the Hedin equation |g^ and investigate the GW- 
approximation |g^ of solid-state physics. 

As a last point, it can be noticed that previous arti- 
cles have presented general planar trees as the structure 
adapted to quantum field theory [p3. In fact, planar trees 
and planar binary trees are equivalent for that purpose. 
Since the number of planar trees with n vertices is equal 
to the number of planar binary trees with 2n — 1 vertices 
0, there is a bijection ^ between planar trees and pla- 
nar binary trees. If i is a planar tree, the bijection can be 
realized recursively by, for example, !?"(•)= • and 



>F{B+{h,t2, ..., tk)) = i^(ti) V >F{B+{t2, 



,tk)), 



where _B_|_ is the grafting operator defined by Kreimer Q, 
0| to build rooted trees. This bijection is also apparent in 
the existence of two methods for the numerical solution of 
differential equations on Lie groups: one based on planar 
trees pq| , the other on planar binary trees |26J . Planar bi- 
nary trees were chosen here because the recursive formulas 
look simpler and because of the mathematical results of 
Loday, Frabetti and collaborators. 
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11 Appendix 

This appendix contains proofs of some of the statements 
contained in the text. 



11.1 Proof of Eq.(||) 

Equation (g) will be proved in two steps. Firstly, it will be 

shown that if 0"(t) satisfies Eq.(|), then ^^^ = '/'^+H0> 
secondly, that the sum over trees is a solution of Eq.(Q). 
The first step is to show that 



'{t;{z}i,n) 
Sv{y) 



^^l+\t;y,{z}i,n). 



(29) 



This will be proved inductively. According to the construc- 
tion of (j)" ( • ) from the initial data A, the relation is true 
for t = • . Now we assume that the relation is true up to 
trees with 2N — 1 vertices. Let i be a tree with 2N + 1 
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vertices. The functional derivative of Eq.(@) gives us 



Sviy) ^ Sv{y) 

(bl-''+\h;z,{z}k+i,n)) + 

n 

^F(</>^-(ii;{^}i,fc), 



At this point intervenes the essential property Eq. m) that 
each tree different from • is generated in a unique way 
by the grafting of two trees ti and ^2- The sum over ti 
and ^2, which are the branches of i, can be replaced by an 
unrestricted sum over all trees ti and t2. Thus 



fe=0 



rn-k+l 



E it2;Z,{z}k+l,n) 

6v{y) 



)• 






A + FiX, 



Since ii and ^2 have less vertices than t, relation ( P9| ) is 
true for them and we obtain 

= 2^F[(Ps {ti;y,{z}iM), 



6X . 
Sv{z)' 



dv{y) 



fc=0 



Cf,]r''+Ht2;z,{z}k+I.n)) + 
n 

^F(/(ti;{z}i,fc), 

fe=0 

0^1;'=+2(i2;2/,^,{z}fe+i,„)),(3O) 

where (/)|j^(i; {z}i_fc) distributes the variables zi and Z2 
over the k positions zi, . . . ,Zk without changing the order 
of Z3, . . . , Zfc. All ways to take two number among k are 
used, so it is clear that (j)^^{t; zi, Z2, ■ ■ ■ , Zk) is symmetric 
in 21,2:2- 

On the other hand, we know from Eq.(g|) that 

r+Ht; y, {2}i^„) = F{cjP{h), 0'i;+2(i2; 2, y, {z},,,)) + 
F(0"+i(ii;y,{2}i,„),0i(i2;2)) + 

n 

J2 F{<l^\h-y, {2}i,fc_i), 0r'^'(^2; 2, {2K„)). 
fc=l 

If we symmetrize y in 0"'*'^, we obtain 
<j,l+\t-y,{z},,^) = F{^\t,),cf>l+i{t2-.z,y,{z}i.n)) + 



A small extension of the previous result is necessary 
to treat the case of QED. The Schwinger equation is now 
X = A + AF{X, j^) and the solution is X = ^^t 0(i), 
where 4>{») is the usual initial data and the recurrence 
relation becomes 



n n—k 
^"^ ^"^ tMl 



fe=Ofe'=0 

F(0^'(ii;{2}fc+i,fc+feO,'/'r'"''+'(^2;2,{z}fe+fc.+i,„)). 

The reasoning is the same as for the simple Schwinger 
equation. We start by proving that -^^7^ = 4>s^^{t). This 
is done by using the recursive definition of cj) to write 
both sides in terms of 0(ti) and (j){t2)- Then we write 
X = '}2n(tf'{t) and we show that, because of the recur- 
rence relation for (/>, it satisfies X = A + AF{X^ Svfz) ^' 

In this paper, the treatment was restricted to linear 
integral operators F. However, the method can treat other 
cases. For instance, if F depends on X^, this square can 
be, in some sense, "linearized" by the planar trees. If A" = 
Et</'°(i),then 



^^ = (E^°w) -EE'^°(^i)<^°(^^) 



tl t2 



A"+l 



iti;y, {z}i.n), (f>Eit2; z)) + 



^.,0 



V-^W, 



E F{cp%ih;y, {z}i,k-i), 0r'+'(^2; 2, {z},,„)) + 

n 

E F{cj,''itu{zh^k), (l>l-i+Ht2:.z, y, {2},+i^„)). 



k=l 



Comparing this with Eq. ( po| ) , we see that the two expres- 
sions are identical, and the property is proved for t. 

If we denote by X the (formal) sum X = J^t'^'^i^)' 
~ '"^ ~ Et'l^Ht-^z). Using 



then from Eq. (|29|) we obtain 
Eq.(H) we can write 

A = E0°w 



Sv{z) 



where ip'^it) (and more generally V'"(i)) is defined by 

n 

r{t) = Y.ci^\hW-\t2). 

k=0 

11.2 Proof of Eq.@ 

Equation ^ will be proved by two methods. Both are 
fast and easy. In the first method, we define a Schwinger 
equation 



b'{') + Y.Fic^"{h)A'it2;z)). 



H 



(31) 



^""^^^ ' 5y N mi r We use property ([52D in the right-hand side : 
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where If this is true up to n, then, from Eq.(n]) 

\ av\^z) dv(z)/ \t\=n+l k=0\ti\=k\t2\=n-k 

This is a Schwinger equation whose solution is given , r- • • /Ctttv 

by Eq.(|, where the map is ^^ deBmtion (Ql, 

IW-Z^ ^if ^^1^+"^ iiV' J' \t\=n+l k=l\t,\=k\t,\=n-k 

which is Eq.(|). XI XI X X*^^!^"*)®^'- 

The upper component of Eg. (pi]) is the first Schwinger k=2 |ti|=fc \t2\=n-k i=i 
equation X = A + F{X, ^^). If we write Y = X + Z 
the lower component is Y — X + G{Y, -^y^)- Therefore 

X + Z is the solution of the composition of equations. The ,^ ,^ ,^ ,^ 

map corresponding to this solution is x(t) = (j}(t) + ii{t). /^ P{t) = ^ /^ 2^ ^^^ ^ • ) ® ^2 + 

The second proof can also be useful. It starts by adding \t\=n+i fc=i |ti|=fc \t2\=n-k 

a parameter s to the Schwinger equations: X{s) = A + n-i n 

sF{X{s)/-^) and Y{s)^X{s) + sG{Y{s)/^). We E E E (^i V 1^.-™) ® l^n- 

take the n-th derivative of Y with respect to s, and we m=ife=m+i |ti|=n-fc 

write F^"^ for its value at s = 0. The chain rule gives _ 

In the second term, we use Eq.([l|) to sum over k by 

k=0 ^ ' ^ ' / ^ J^-/c V Yk-m = J^-m+1 — J^-m V • . 

k—rn+1 

It can be shown recursively that 

Therefore 

l*l=" ^ P{t) = T. E E (^iV.)®t2 + 

where ip satisfies Eq.(||). The result follows by expanding |t|=n+i fc=i |ti|=fc \t2\=n-k 

X{s) and Y{s) as power series of s and taking its value n-i 

for s = 1. E E E *i ® *2 - 

m=l |ti|=n— m+1 \t2\=m 

11.3 Proof Of Eq.® g ^ (iiV.)®i2 

Let X(5) = l + y]^^* 5l*l(/)(t). Let US show that the series "^~ lUI-n-m 

for Y{s) = 1 - 1/X(s) is given by Y{s) = ^m . «'*'^(*)' =5^ J] ^ ti ® ^2- 

where V'lO is defined by Eq.(|2^). m=i jtj|^„+i_™ i^^i^^ 

The first step is to prove that, if P{t) is defined by 

Eg. (psi) , then, with an abuse of notation. To complete the proof, we start from X{s) = 1 + 

^j / , sl*l0(i) and we define Y{s) — J2t=^ • ^'''■'/'(O; where 

P(r„) ^y^Yk^Yn-k (32) ^W is given by V'l • ) =-0, ipit) = 0(t) - ^ 0(Mi) «) V(w*)- 

^ Thus 

or more precisely _^ ,,, _^ ,,, _^ i,, ^ 

E *"^w = E *"<^w - E ^" E <^(^^)^("')- 
^w = E(E")®( E -)■ ^ 

|t|=n fe=i |«|=fc |i;|=ri-fc From Eq. (|32[) , we see that |i| = |ui| + \vi\ and 



As usually, this will be proved recursively. The property v^ \t\ , ^ ,\ _ v^ \*\A,(f') _ 

is true for n = 2, because ■^ ^ 

\t\>\ \t\>i 

V V l"l>0 |f|>0 
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From the definition of X{s) and Y{s) we deduce 

Y{s) ~ ^(V) - X{s) - 0(V) - 1 - iX{s) - l)Y{s). 

From the definition of ^{t) we get V'(\^) = 0(\^)7 and 
Y{s) satisfies the equation 

Xis)-Xis)Y{s)^l, 

or X{s) — 1/(1 — Y{s)). To simplify the notation, we have 
assumed (/>(•) = 1. The general case is proved similarly 
and leads to Eq.(|2^). 
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